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ABSTRACT 

Several exoplanets have recently been imaged at wide separations of >10 AU from their parent 
stars. These span a limited range of ages (<50 Myr) and atmospheric properties, with temperatures 
of 800-1800 K and very red colors (J — H > 0.5 mag), implying thick cloud covers. Furthermore, 
substantial model uncertainties exist at these young ages due to the unknown initial conditions at 
formation, which can lead to an order of magnitude of uncertainty in the modeled planet mass. Here, 
we report the direct imaging discovery of a Jovian exoplanet around the Sun-like star GJ 504, detected 
as part of the SEEDS survey. The system is older than all other known directly-imaged planets; as 
a result, its estimated mass remains in the planetary regime independent of uncertainties related to 
choices of initial conditions in the exoplanet modeling. Using the most common exoplanet cooling 
model, and given the system age of 160lgo° Myr, GJ 504 b has an estimated mass of 4 ^'q Jupiter 
masses, among the lowest of directly imaged planets. Its projected separation of 43.5 AU exceeds 
the typical outer boundary of ^30 AU predicted for the core accretion mechanism. GJ 504 b is also 
significantly cooler (510^20 K) and has a bluer color ( J — H = —0.23 mag) than previously imaged 
exoplanets, suggesting a largely cloud-free atmosphere accessible to spectroscopic characterization. 

Thus, it has the potential of providing novel insights into the origins of giant planets, as well as their 
atmospheric properties. 

Subject headings: planetary systems — stars: individual (GJ 504) — planetary systems: formation 
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1. INTRODUCTION 

More than 890 extrasolar planets (or exoplanets) have 
been discovered by now (Schneider et al. 2011), and 
NASA’s Kepler mission alone has recently added more 
than 2,000 exoplanet candidates (Batalha et al. 2013). 
The vast majority of exoplanets have been discovered 
using indirect detection techniques that infer the pres- 
ence of the exoplanet by monitoring the host star, such as 
planetary transits or radial- velocity variations. The tech- 
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niques require observations over at least one orbital pe- 
riod, making them impractical for detecting long-period 
exoplanets at large separations. Meanwhile, in spite of 
the technical challenges, 8-m class telescopes equipped 
with adaptive optics systems and/or coronagraphs have 
recently reached the high-contrast performance neces- 
sary to image massive planets at large orbital separations 
(Marois et al. 2008, 2010; Lagrange et al. 2010). Direct 
imaging surveys using such instruments provide measure- 
ments of an exoplanet’s position and luminosity, which 
can be used to estimate its mass using the age of its host 
star and models of luminosity evolution for planetary- 
mass objects. 

Because massive exoplanets cool and fade with 
time, direct-imaging searches have been most successful 
around young stars. Mass estimates for young planets de- 
pend strongly on the assumed system age, planetary at- 
mosphere models, and initial thermodynamic state, lead- 
ing to large uncertainties in the inferred mass. Plane- 
tary mass estimates conventionally assume a “hot start” , 
in which the planet is initially in a high-temperature, 
high-entropy (and hence, luminous) state (Baraffe et al. 
2003; Burrows et al. 1997). However, recent theoret- 
ical models suggest that giant planets produced ac- 
cording to standard formation theories could initially 
be much colder (Mar ley et al. 2007; Fortney et al. 2008; 
Spiegel & Burrows 2012). The difference between the 
mass-luminosity relations of the “hot-start” and “cold- 
start” models decreases with exoplanet age; the models 
converge after ^100 Myr for a 5 Mj up planet, and ~1 Gyr 
for a 10 Mj up planet (Marley et al. 2007; Fortney et al. 
2008; Spiegel & Burrows 2012). 

Previously imaged exoplanets are all younger than 
50 Myr (Marois et al. 2008; Lagrange et al. 2010; 
Carson et al. 2013). The hot-start models imply plan- 
etary masses (Marois et al. 2008, 2010; Lagrange et al. 
2010), while cold-start models cannot reproduce 
the high observed luminosities (Marley et al. 2007; 
Spiegel & Burrows 2012). Therefore, their planetary sta- 
tus would be negated if the cold-start models apply at 
such young ages. Nevertheless, independent constraints, 
such as dynamical stability analysis in the case of the 
HR 8799 system (e.g., Fabrycky & Murray-Clay 2010; 
Marois et al. 2010; Sudol & Haghighipour 2012), can in 
some cases be used to exclude the coldest range of ini- 
tial conditions, and confine their masses to the planetary 
mass regime. In addition, we note that neither hot-start 
nor cold-start models have been calibrated against dy- 
namical mass estimates. These arguments suggest that 
no fully established model for estimating the masses of 
directly imaged planets exists as of yet. However, direct 
imaging of an older system with a well-determined age 
could remove much of the model uncertainty related to 
the unknown initial conditions. 

A common feature of the previously imaged planets 
is large effective temperatures (> 800 K; Marois et al. 
2008), and very red colors implying the presence of thick 
clouds in their atmospheres. Furthermore, they have 
generally been discovered as companions to massive (> 
1.5 Msun) stars (e.g., Baines et al. 2012; Lagrange et al. 
2010; Carson et al. 2013), hence they populate a limited 
range in atmospheric conditions and host star properties. 
Widening this range would greatly enhance our under- 
standing of the population of exoplanets in wide orbits. 


Here, we report the detection of an exoplanet around 
the nearby Sun-like star GJ 504 using high-contrast near- 
infrared imaging on the Subaru Telescope. The projected 
separation of the planet (GJ 504 b) is measured to be 
43.5 AU from the star, and the mass of planet (GJ 504 
b) is estimated to be a few Mj up ; because the system has 
an age of 160lgQ° Myr, its mass estimate is only weakly 
dependent on the uncertainty of initial conditions in the 
cooling models of giant planets. In Section 2, we describe 
the properties of GJ 504. In section 3, the observations 
and data reductions are detailed. Section 4 presents a 
derivation of the mass of GJ 504 b and other interest- 
ing properties of this planet. Section 5 investigates the 
atmospheric properties and possible origin of GJ 504 b. 
Our results are summarized in Section 6. 

2. GJ 504 STELLAR PROPERTIES 

GJ 504 is an isolated GO-type main-sequence star. It 
has a mass of 1.2 Ms un and an age of 160lgQ° Myr. In 
the following sections, we will describe the stellar prop- 
erties in more detail, including its age estimate. Table 1 
summarizes the properties of GJ 504. 

2.1. Age of GJ 504 

For an isolated main-sequence star like GJ 504, sev- 
eral methods of age estimation can be imagined, includ- 
ing gyro chronology, chromospheric or coronal activity, 
lithium abundance, kinematics, or isochrones. Among 
them, two particularly relevant age indicators for GJ 504 
are gyrochronology based on stellar rotation, and chro- 
mospheric activity as traced by Ca II H and K emission, 
and this study thus use these indicators as the fiducial 
age estimators. Both stellar rotation and activity de- 
cline as a Sun-like star ages and sheds angular momen- 
tum in its stellar wind, and both age indicators have 
been accurately calibrated using stars in open clusters 
(Mamajek & Hillenbrand 2008, hereafter MH08). Chro- 
mospheric activity is powered by the stellar magnetic 
dynamo, which is closely related to the stellar rotation 
rate. Gyrochronology, which directly measures this rate, 
is therefore the most direct, and may be the most reli- 
able, of our age estimators (MH08), and allows us to ob- 
tain a best estimate of 160lgQ Myr for the age of GJ 504. 
However, adopting the age as given by the star’s chromo- 
spheric activity in addition to that from gyrochronology, 
we conservatively adopt 160lgQ° Myr as the age estimate 
of GJ 504 system (see Fig. 1 for summary). We describe 
these age estimations in detail below, along with several 
other age indicators that are less well suited for the pur- 
pose. 

2.1.1. Gyrochronological Age 

The rotation of a star slows as the star ages and 
the stellar wind carries away angular momentum. Gy- 
rochronology, which estimates a star’s age from its rota- 
tion rate and B — V color (i.e., mass), is generally consid- 
ered to be the most direct, and possibly the most accu- 
rate, method for inferring stellar age (MH08). Chromo- 
spheric and X-ray activity are powered by the stellar dy- 
namo, which is closely related to the star’s rotation rate. 
Indeed, chromospheric age estimators often use stellar 
activity to infer the rotation period, which is then used to 
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TABLE 1 

Properties of GJ 504 system 


Property 

Primary 

Planet 

Spectral type 

G0V a 

late T - early Y 

Distance (pc) 

17.56 ± 0.08 b 

Proper motion (mas/yr); up 

-338.83 ± 0.25, 

190.24 ± 0.21 b 

Effective temperature; T e ff (K) 

6234 ± 25 c 

510l“ d 

Surface gravity; logg 

4.60 ± 0.02 c 

q q+0.4d 

o.n_ 0 2 

Luminosity; log (L/L®) 

0.332 ± 0.032 c 

-6 09 + °'° 6d 

o.un_ o os 

Metallicity; [Fe/H] 

0.28 ± 0.03 c 

— 

Adopted age (Myr) 

i6o; 

T350 

-60 

Age by gyrochronology (Myr) 

160 

1 + 
C^> -4 

o o 

Age by chromospheric activity (Myr) 

330 ± 180 

Mass 

1.22 ± 0.08 M© 

4.0lto M Ju P d 

J-band mass 

— 

3.5i? ® M Jup 

iL-band mass 

— 

3-Oto's M Ju P 

ATg-band mass 

— 

6.5±®’° M Jup 

I/-band mass 

— 

3.5±?'° M Jup 

NIR brightness 6 (mag) (J; ~1.2 /xm) 

4.106 

19.91 ± 0.15 

(if; ~1.6 jiim) 

3.860 

20.14 ± 0.14 f 

(K s \ ~2.2 /im) 

3.808 

19.19 ± 0.28 

(L x ; ~3.8 jum) 

3.94 

16.84 ± 0.19 


a Valdes et al. (2004), Takeda et al. (2005) 
b (van Leeuwen 2007) 
c Valenti & Fischer (2005) 

d GJ 504 b’s effective temperature, surface gravity, bolometric luminos- 
ity, and mass are estimated by comparing its age (160igo° Myr) and 
JHK S L ' -band magnitude with the most commonly used cooling model 
(Baraffe et al. 2003). For quantifying best estimates and lower/upper 
limits for all of these quantities, we apply the same calculation as de- 
scribed in section 4.2. 

e The provided errors of the stellar JHK photometry are less than 0.01 
mag (Kidger V Martm-Luis 2003), while L' photometry was obtained 
in this work with an error of 0.09 mag (see section 3.5). 
f The if -band photometric magnitude is based on the weighted mean of 
the March and May values. 


infer an age (MH08). The rotation period of GJ 504 has 
been directly measured to be 3.328 days (Messina et al. 
2003) and 3.33 days (Donahue et al. 1996). We adopt the 
average of these measurements, 3.329 days, as the star’s 
rotation period. The possible error range of < 0.01 days 
is small enough not to affect the uncertainty of the age 
estimate. 

The gyrochronology can be applied to I-sequence stars 
whose dynamo mechanism is presumed to originate at 
the interface between the convective and radiative zones 
in the stellar interior (Barnes 2007). In contrast, this 
method cannot be applied to C-sequence stars that are 
fully convective and may thus have a different dynamo 
mechanism that prevents efficient spin-down. We can 
judge that GJ 504 is an I-sequence star using two dif- 
ferent sets of criteria specified in the literature (Barnes 
2007; Dupuy et al. 2009): one in which we relate the age- 
normalized rotation of GJ 504 to other stars of similar 
B — V colors, and the other in which the X-ray luminos- 
ity and Rossby number are compared to other chromo- 
spherically active stars. Both of these methods confirm 
that GJ 504 is consistent with an I-sequence star. This 
conclusion agrees with the previous classification (Barnes 
2007) of GJ 504 as an I-sequence star. 

To rigorously establish the age of GJ 504, we use 
three gyrochronology relations independently proposed 


(MH08; Barnes 2007; Meibom et al. 2009). All the re- 
lations determine a stellar age as a function of rotation 
period and B — V color. Using the rotation period (3.329 
days) and B — V color (0.585 db 0.007) of GJ 504, we de- 
rive ages of 140-180 (MH08), 100-140 (Barnes 2007), and 
110-230 Myr (Meibom et al. 2009), respectively. The er- 
ror ranges for the age estimates were derived based on 
all the possible sources of error: the observed scatter 
in the period-age relation (0.05 dex; MH08) or errors 
of coefficients in the period- age function (Barnes 2007; 
Meibom et al. 2009), and the error of the B — V color. 
Conservatively adopting all these ranges 1 , we estimate 
the age of GJ 504 to be 160tgo Myr. 

In order to verify the adequacy of GK 504’s age 
assessed from the empirical relations of gyrochronol- 
ogy (see above), we directly compare the rotation pe- 
riod and B — V color of GJ 504 to those of stars in 
young clusters (see Fig. 2). The rotation data for 
stellar members in Pleiades (130 Myr; Hartman et al. 
2010), M35 (150 Myr; Meibom et al. 2009), M34 (220 

1 For the optimal value of the gyrochronological age estimation, 
we adopt 160 Myr which was calculated by MH08, because MH08 
have updated the gyrochronological relation of Barnes (2007). 
Both relations of MH08 and Barnes (2007) have been derived by 
utilizing several stellar clusters, while Meibom et al. (2009) used 
only a single cluster. 
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Fig. 1. — GJ 504’s ages independently inferred from gyrochronol- 
ogy, chromospheric and X-ray coronal activity, or lithium technique 
are summarized together with our adopted age for the system. 
The squares indicate the optimum estimates based on each tech- 
nique, and the horizontal bars correspond to their error ranges. Be- 
cause the gyrochronology and chromospheric activity could be the 
most reliable age indicators for isolated G-type stars with evolution 
stages near ZAMS, we use those as the fiducial age estimators. The 
best-estimate in our adopted age for GJ 504 is the gyrochronology 
technique, which directly traces the stellar rotation than the activ- 
ity. The estimates from lithium or X-ray activity are not adopted 
but consistent to those from others. 

Myr; Meibom et al. 2011), Mil (230 Myr; Messina et al. 
2010), M37 (550 Myr; Hartman et al. 2009), Hyades (630 
Myr; Radick et al. 1987, 1995; Delorme et al. 2011), and 
Praesepe (600 Myr; Delorme et al. 2011) can be used for 
this purpose. As shown in Fig. 2, the stars with simi- 
lar B — V colors to GJ 504 in 130-230 Myr old clusters 
have periods that are the most consistent with that of GJ 
504, and the rotation periods in the older 500-600 Myr 
cluster are significantly longer 2 , thus verifying the age 
range of GJ 504 derived above. Additionally, combin- 
ing the rotation data of stars in Mil with those in M34, 
Messina et al. (2010) derive a median rotation period of 
G-type stars in these clusters of 4.8 days, again support- 
ing the age estimate of 160^60 Myr for the G-type star 
GJ 504 (3.329 day period). 

Since gyrochronology is a largely empirical method, 
there could in principle be uncertainties involved in the 
methodology that are not accounted for in the calibra- 
tions. For instance, if we hypothesize that the calibra- 
tion could be metallicity-dependent, this would lead to 
systematic errors for targets with significantly different 
metallicity from the calibration clusters. For the case 
of GJ 504, it has a non-extreme value in this regard, 
hence it is unlikely to be subject to such uncertainties. 
Another possible systematic uncertainty in gyrochronol- 
ogy involves the fact that rotation is easier to detect in 
rapidly rotating systems, which could result in a bias in 

2 In addition to I-sequence stars, Fig. 2 includes rotation data of 
C-sequence stars, which are widely scattered in the plot. However, 
the scatter is much smaller for the older clusters. GJ 504 is best 
represented by I-sequence stars, and young I-sequence stars such 
as provided in M35 or M34 (Meibom et al. 2009, 2011) are indeed 
in good agreement with GJ 504. 



Fig. 2. — Stellar rotation of GJ 504, and of stars with B — V < 1.0 
in various clusters. The horizontal axis shows the B — V colors of 
the stellar samples, while the vertical axis denotes their rotations. 
The data of stars in ~100-200 Myr old clusters (Pleiades, M35, 
M34, and Mil) are plotted as open circles, and those of ~500-600 
Myr old clusters (M37, Hyades, and Praesepe) are plotted as blue 
squares. See the main text for the references of these data. The un- 
certainties in stellar rotation periods should be smaller than ~10% 
in general (cf. Donahue et al. 1996; Barnes 2007; Hartman et al. 
2009). The squares with B — V colors bluer than 0.55 mag are 
M37 stars, with uncertainties of ~0.01 mag. Thus, GJ 504 does 
not overlap with stars in 500-600 Myr clusters. 

the empirical gyrochronology relation. If this were a sig- 
nificant effect, it would cause the apparent mean rotation 
in a cluster of a given age to appear more rapid than it is 
in reality. Hence, age estimates for individual stars based 
on a relation calibrated against such clusters would re- 
sult in an over-estimation (i.e., an individual star would 
appear older than it really is). Given that no other age 
indicators imply a significantly younger age for GJ 504 
this is also unlikely to be a real issue, but nonetheless, 
such uncertainties are important to keep in mind in the 
age analysis. 

2.1.2. Chromospheric Activity Age 

Stellar chromospheric activity, as traced by Ca II H 
and K emission, is a good indicator of stellar age. The 
strength of Ca II H and K emission is parametrized by 
R! l ik, roughly the power in the H and K lines normal- 
ized to that in the underlying photospheric continuum. 
A relation between R'h K , B — V color, and age has re- 
cently been derived for solar- type stars (MH08). For GJ 
504, we adopt the value of log R'h K = — 4.45 measured 
at Mt. Wilson Observatory with a temporal baseline of 
30 years (Radick et al. 1998), the longest such baseline 
in the literature. In order to minimize the uncertainty in 
R'h K , the value for which the star has been monitored for 
the longest possible time should be adopted (see MH08), 
which also protects against accidentally adopting an out- 
lier in the activity variation. 

This activity level may be used to infer a Rossby 
number, the ratio of the stellar rotation period to the 
timescale of convective overturn, of 0.60 =b 0.10 following 
MH08 3 . We use this result 4 , together with a B — V color 
of 0.585 =b 0.007 (van Leeuwen 2007), to estimate an age 
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of 330 ±100 Myr. We further add a 0.2 dex scatter in 
this relation, estimated from clusters 5 (MH08). We as- 
sume this scatter to be uncorrelated with other sources 
of error. As a result, we obtain a chromospheric age 
estimate of 330 ± 180 Myr. 

Other studies using data from Mt. Wilson Observa- 
tory have measured R'h K values of -4.443 (Baliunas et al. 
1996) and -4.486 (Lockwood et al. 2007), which differ 
slightly from our adopted value of -4.45. The chromo- 
spheric activity level was also measured to be R'u K = 
—4.44 using 14 years of observations at Fairborn and 
Lowell Observatories (Hall et al. 2009). Many of these 
studies’ observational baselines overlap with the source 
of our adopted value (Radick et al. 1998). Adopt- 
ing the mean of these measurements would have little 
effect on our chromospheric age of 330 ±180 Myr. 

A primary uncertainty in age determination from chro- 
mospheric activity is long-term activity cycles, which 
may easily exceed the timescales over which the activity 
is measured, and thus may potentially provide a system- 
atic error if the activity is measured at an extreme but is 
taken to represent the mean activity for the star. In addi- 
tion, as in the case of gyro chronology, the stellar activity 
also depends on the metallicity, possibly contributing to 
the uncertainty. Since activity as an age estimator is 
thought to be closely related to gyrochronology through 
the stellar dynamo, the possibility of measuring both ro- 
tation (which is robust against activity cycles) and activ- 
ity in the same star as for the case of GJ 504 can alleviate 
this uncertainty. 

2.1.3. X-ray Activity 

Young solar- type stars are active X-ray emitters. As 
for Ca II H and K, a relation has been derived (MH08) 
between stellar age, color, and X-ray luminosity. The 
best estimate of the X-ray luminosity of GJ 504 is from 
the ROSAT all-sky survey of nearby stars (Hiinsch et al. 
1999), which reports a ratio of X-ray to bolometric power 
logLx/Tboi = —4.42. This X-ray activity may be used 
to infer a Rossby number of 0.54 ± 0.25 (MH08) and 
lead to the age estimate of 90-530 Myr for GJ 504. Be- 
cause the temporal baseline for the X-ray measurement 
is much shorter than for the multi-decade measurements 
of Rhk? inferred Rossby number is much more un- 
certain. Therefore, we do not adopt this estimate based 
on the X-ray emission, although it is consistent with the 
estimate from other methods. 

2.1.4. Lithium absorption 

The strength of lithium absorption lines declines as 
a star burns its primordial supply of lithium, provid- 
ing another indicator of stellar youth (Soderblom 2010). 
However, many effects, such as variations in the pri- 
mordial supply of lithium, the extreme sensitivity of 

measurement of to the Rossby number. The error in the in- 

ferred Rossby number was measured to be 0.10 for the multi-decade 
Mount Wilson data (see section 4.1.2 and summary of MH08). 

4 The Rossby number is converted into a rotation period, which 
is used in the gyrochronology relation (i.e., activity-rotation-age 
relation) to estimate the stellar age. This scheme improves the ac- 
curacy of age estimation from chromospheric activity (see MH08). 

5 MH08 infer the uncertainty of their activity-rotation-age rela- 
tion by calculating the rms of residuals between the ages of clusters 
and those of their stellar members, which were estimated by ap- 
plying the activity-rotation-age relation to the members. 


early lithium burning to the stellar accretion history 
(Baraffe & Chabrier 2010), and the importance of non- 
local thermodynamic equilibrium abundance corrections 
(Takeda & Kawanomoto 2005), make it difficult to ob- 
tain a reliable age estimate from lithium abundances. 
Therefore, we do not employ lithium as an indicator to 
date GJ 504, but it can provide a clue to the age of GJ 
504 as discussed below. The equivalent width of Li in G J 
504 has been measured to be 0.08 A, implying a lithium 
abundance of log n(Li) = log(N(Li)/N(H)) ± 12 = 2.9 
(Takeda & Kawanomoto 2005). Extensive observations 
of well-dated open clusters (Sestito & Randich 2005) en- 
able a lithium age estimate for GJ 504. At the effective 
temperature of a GO star (^6000 K), these cluster esti- 
mates can be used to infer an age range of approximately 
30-500 Myr, fully consistent with our gyrochronological 
and chromospheric age estimates. 

2.1.5. Kinematics 

GJ 504 is not currently identified as a member of any 
clusters or moving groups, making it impossible to use 
population-based age estimators. The kinematic velocity 
of GJ 504 (U=- 38.6, V=1.5, W= 18.0; Mishenina et al. 
2004) is consistent with thin disk stars, thus providing 
no tight constraint on the age, but remaining consistent 
with the young age provided by gyrochronology and chro- 
mospheric activity. 

2.1.6. Isochrones 

Since stars evolve along mass tracks in an HR-diagram 
with age, a common age determination method is based 
on isochronal analysis. This method can provide rea- 
sonable estimates for old stars that have evolved signif- 
icantly off the main sequence and very young stars that 
are still undergoing rapid contraction. However, it is 
generally not useful for stars near the zero-age main se- 
quence, around which stellar evolution proceeds slowly 
(e.g., Soderblom 2010). Indeed, generic isochronal esti- 
mates in the literature are known to give erroneous ages 
for young stars by up to two orders of magnitudes. 

As illustrative examples, we consider the cases of 
HD 152555, HIP 10679, and HD 206860. These are all 
early G-type (G0-G2) stars which are members of young 
moving groups (MGs), hence they are similar types of 
stars to GJ 504 with well-determined young ages, and 
thus can be used for comparison to isochronal ages of 
e.g., Valenti & Fischer (2005). HD 152555 is a member 
of the AB Dor MG (e.g., Torres et al. 2008) with an age 
of 50-100 Myr, whereas Valenti & Fischer (2005) give 
an age of 4.0 Gyr. HIP 10679 is a member of the (3 Pic 
MG (e.g., Torres et al. 2008) with an age of ^10 Myr, 
whereas Valenti & Fischer (2005) give an age of 3.0 Gyr. 
Finally, HD 206860 is a member of the Her-Lya MG (e.g., 
Lopez- Santiago et al. 2006) with an age of ^200 Myr, 
whereas Valenti & Fischer (2005) give an age of 3.1 Gyr. 

In light of these large discrepancies, it is unsurpris- 
ing that literature estimates for the isochronal age of 
GJ 504 vary strongly from hundreds of Myrs to several 
Gyrs (e.g., Valenti & Fischer 2005; Takeda et al. 2007; 
Holmberg et al. 2009; da Silva et al. 2012). This implies 
that isochronal analysis should simply be disregarded 
for a star such as GJ 504, compared to the more reli- 
able age estimators listed above. Nonetheless, since some 
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Fig. 3. — Isochronal analysis of GJ 504. Left: Reproduction of the da Silva et al. (2012) result. Also plotted is another data point where 
all values are kept equal, except for the temperature which has been replaced with the one in Valenti & Fischer (2005), illustrating the 
real uncertainty in temperature. With the Valenti & Fischer (2005) temperature, the data is fully consistent with the ages derived from 
gyrochronology and activity. Right: The same analysis, again with all quantities except for temperature fixed to the da Silva et al. (2012) 
values, but replacing luminosity with surface gravity. In this framework, >3 Gyr ages are disfavored regardless of which temperature is 
chosen. 


of the isochronal literature estimates provide Gyr-range 
ages that are seemingly inconsistent with the younger 
age range we have adopted for GJ 504, we discuss their 
application to this star in more detail below. 

An important limiting factor in isochrone analysis for 
age determination of an isolated main-sequence star such 
as GJ 504 is the sensitivity to uncertainties in input stel- 
lar parameters. This is in particular true for properties 
such as effective temperature and surface gravity, where 
the uncertainties are model-dependent and often tend 
to be under-estimated in the literature. As an illustra- 
tion, we consider the isochronal analysis in da Silva et al. 
(2012) based on the Yonsei-Yale isochrones (Yi et al. 
2001), which gives an age of ^3 Gyr, and appears to 
be inconsistent with ages of <500 Myr within their 
quoted error bars. We have reproduced the analysis of 
da Silva et al. (2012) and confirmed these results with 
the same methodology and input observables. This is 
plotted in Fig. 3. Also plotted in the same figure is 
the exact same analysis, but using the temperature esti- 
mates of GJ 504 from Valenti & Fischer (2005) instead; 
note that all other values are kept as in da Silva et al. 
(2012). This comparison immediately illuminates one 
of the considerable uncertainties related to isochronal 
analysis, which is uncertainties in the temperature scal- 
ing. The two temperature estimates in Valenti & Fischer 
(2005) and da Silva et al. (2012) are inconsistent within 
their quoted error bars, implying that neither error es- 
timate should be trusted when translating this into an 
error on the age from the isochronal comparison. More- 
over, we find that while the temperature estimate of 
da Silva et al. (2012) gives the aforementioned age of 
^3 Gyr and is inconsistent with <500 Myr ages, the 
temperature estimate of Valenti & Fischer (2005) is fully 
consistent with these younger ages, and inconsistent with 
the rsj 3 Gyr age. 

It is also important to note that due to pre-main 
sequence evolution, there is some overlap between the 
isochrones of stars arriving at the main sequence, and 
those leaving it. For instance, the 100 Myr isochrone 
overlaps with the 1 Gyr isochrone in the relevant temper- 


ature range. Hence, while the age is quoted as ~1.4 Gyr 
in Valenti & Fischer (2005) which is reproduced in Fig. 
3, it can be equally well interpreted as being consistent 
with our adopted age range, and in fact preferentially 
supporting the youngest ages. 

Furthermore, while the temperature and luminosity in 
da Silva et al. (2012) gives a certain age from compar- 
ison with the Yonsei-Yale isochrones, this does not ac- 
count for the relatively high surface gravity they quote. 
We can see this in the second panel of Fig. 3, where 
again we perform an isochronal analysis, but in surface 
gravity versus temperature as opposed to luminosity ver- 
sus temperature. Here, regardless of which temperature 
value used, the data are consistent with a <500 Myr age, 
but inconsistent with a ^3 Gyr age. 

It is instructive to think of these discrepancies in terms 
of the radius of the star. If we adopt the da Silva et al. 
(2012) temperature, then the luminosity is high for its 
mass, implying a large radius. On the other hand, 
the surface gravity is also high for its mass, implying 
a small radius. Hence, there is a tension in the ra- 
dius between the luminosity and surface gravity if the 
da Silva et al. (2012) temperature is adopted, such that 
no self-consistent picture can be reached by either in- 
creasing or decreasing the stellar radius. By contrast, 
if we adopt the Valenti & Fischer (2005) temperature, 
then this tension disappears, and a self-consistent picture 
can be reached between the temperature, luminosity, and 
surface gravity - and simultaneously, we get a predicted 
age which is consistent with the estimates given by the 
age estimates in the previous sections. This could there- 
fore be interpreted as supporting the higher temperature 
of Valenti & Fischer (2005) and thus a younger age. 

Hence, we note that the isochronal age is readily con- 
sistent with gyrochronology and activity, using values 
for temperature, luminosity, and surface gravity pub- 
lished in the literature. This is also consistent with the 
arguably most detailed Bayesian isochronal analysis of 
Takeda et al. (2007), which was a study meant to im- 
prove on the analysis of Valenti & Fischer (2005) using 
the same input data, and which gives an age limit of 
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<760 Myr. Nonetheless, the above discussions imply 
that it is not feasible to provide age constraints from 
isochronal dating that are as meaningful nor as stringent 
as those of gyrochronology and chromospheric activity. 
Thus, we do not include it in our adopted age range, but 
merely conclude that it is broadly consistent with such 
an age. 

Finally, it can be noted that the existence of a star with 
such a short rotation period as GJ 504 (~3 days) would 
be unprecedented at the Gyr-scale ages predicted by 
some isochrone analyses, such as da Silva et al. (2012). 
Indeed, although Reiners & Giampapa (2009) found an 
unusually rapid rotator in the ~4 Gyr old M67 cluster, 
it has a rotation of approximately half the Solar rota- 
tion rate (P© ~ 30 days), which is significantly slower 
than the rotation of GJ 504. Aside from known close 
spectroscopic binaries, it is the largest known outlier at 
such an age. No indications of chromospheric activity 
corresponding to such rapid rotation have been found 
among 60 stars in the M67 cluster (Giampapa et al. 
2006), which supports the notion that rapid rotators like 
GJ 504 must be exceedingly rare at ages of a few Gyr. 

2.1.7. Summary of Age Estimate for GJ 504 

As discussed above, we adopt the age estimated with 
gyrochronology as the most likely age ( 160+60 Myr) for 
GJ 504. However, we conservatively adopt an age of 
160 ±lg° Myr, with chromospheric activity providing the 
upper bound and gyrochronology providing the lower 
bound. Age indicators based on the Lithium abundance, 
X-ray luminosity, kinematics, and isochrone technique 
are less precise than gyrochronology and chromospheric 
activity: in addition, for the Li-age estimate it is more 
difficult to quantify the uncertainty. We therefore do 
not attempt to combine those with the age range esti- 
mated from the gyrochronological and chromospheric ac- 
tivity techniques, although their estimates are consistent 
with our adopted age of GJ 504. Future asteroseismology 
studies, or detailed kinematic studies identifying GJ 504 
as a member of a moving group, may be able to refine 
our age estimate. The estimated ages are summarized in 
Table 1 and Fig. 1. 

2.2. Other Information 

GJ 504 was previously reported to have a stellar 
companion at a projected separation of about 650 
AU (Poveda et al. 1994), but subsequent observations 
showed that it does not share a common proper motion; 
the GJ 504 proper motions are -333.83 d= 0.25 mas/yr 
and 190.24 db 0.17 mas/yr in right ascension and decli- 
nation (van Leeuwen 2007), respectively, while those of 
the stellar companion candidates are 22.9 db 11.2 mas/yr 
and -51.7 db 11.2 mas/yr in each direction (Roeser et al. 
2010). Hence, we identify the stellar companion candi- 
date as an unrelated background star. 

The effective temperature, surface gravity, and metal- 
licity [Fe/H] of GJ 504 was inferred by directly fit- 
ting the observed spectrum to the synthetic spectrum 
(Valenti & Fischer 2005). Interestingly, GJ 504 has a 
super-solar metallicity [Fe/H] =0.28 (Valenti & Fischer 
2005). In order to estimate the mass of the central star, 
we use an empirical relation of Torres et al. (2010), which 
is a function of the surface gravity, effective temperature, 


and metallicity of GJ 504 (Valenti & Fischer 2005). The 
error of the estimate was provided by including uncer- 
tainties in the spectroscopically determined input pa- 
rameters, the reported intrinsic scatter in the relation, 
and correlations of the best-fit coefficients (Fleming et al. 
2010). As a result, the mass of GJ 504 is estimated to 
be 1.22 dz 0.08 solar masses. The mass of GJ 504 can 
also be derived based on matching its spectroscopically 
determined surface parameters with theoretical stellar 
evolution models (Takeda et al. 2007; Valenti & Fischer 
2005). The theoretical models infer the mass of GJ 504 
to be 1.29io;°4 (Takeda et al. 2007) or 1.28 ± 0.03 M Sun 
(Valenti & Fischer 2005), overlapping with the empirical 
estimate. For our analyses, we employ the mass estimate 
based on the empirical function for GJ 504. 

Broad band observations of the infrared excess above 
the photosphere at 24 microns and 70 microns place con- 
straints on the distribution of dust in the inner (< 10 AU) 
and outer disks. A recent Spitzer survey showed that 
4% of Sun-like stars have 24 micron excess, while 16% 
of Sun-like stars have 70 micron excess (Trilling et al. 
2008). GJ 504 was shown to not have excesses at ei- 
ther of these band passes (Tdust / ^star < 2.1 x 10 -6 ), and 
therefore does not seem to harbor a significant debris 
disk (Bryden et al. 2006). 

3. OBSERVATIONS AND DATA REDUCTION 

GJ 504 was observed as part of the Strategic Explo- 
ration of Exoplanets and Disks with Subaru (SEEDS) 
survey (Tamura 2009), which aims to detect and char- 
acterize giant planets and circumstellar disks using the 
8.2-m Subaru Telescope. In 2011 and 2012, we ob- 
tained J- (rsj 1.2 /im), H- (^1.6 /im), and K s - (~2.1 
/im) band images using the High Contrast Instrument 
for the Subaru Next Generation Adaptive Optics (Hi- 
CIAO; Suzuki et al. 2010) with A0188, a 188 actua- 
tor adaptive-optics system (Hayano et al. 2008). Fur- 
thermore, we obtained follow-up L'-band (^3.8 /im) im- 
ages using InfraRed Camera and Spectrograph (IRCS; 
Kobayashi et al. 2000) and A0188. See Table 2 for the 
summary of GJ 504 observations. 

3.1. HiCIAO Observations 

We discovered GJ 504 b in iL-band observations using 
HiCIAO on March 26, 2011 with A0188. We used GJ 
504, its host, as the natural guide star, and used an at- 
mospheric dispersion corrector (ADC) to prevent it from 
drifting on the detector due to differential refraction in 
the visible and near-infrared. We used a circular occult- 
ing mask 0'/4 in diameter and angular differential imag- 
ing (ADI; Marois et al. 2006) to remove the starlight and 
the stellar speckles. ADI keeps the telescope pupil fixed 
relative to the camera, while the field-of-view rotates as 
the target moves across the sky. As a result, diffraction 
and speckle patterns in the stellar point spread function 
(PSF), caused by the major part of telescope optics, re- 
main fixed on the detector while real on-sky sources, such 
as a planet, rotate about the natural guide star. This al- 
lows ADI data processing to distinguish between stellar 
speckles and real objects, and to remove the starlight 
while preserving potential planets. 

We conducted several follow-up observations in 2011 
and 2012 using the same instrument and configuration. 
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TABLE 2 

Observation Log and Astrometric Measurements of GJ 504 b. 


Obs. Date 
(UT) 

Inst. 

Band 

Total DIT 
(min) 

Mode 

FRA 

(°) 

Proj. Sep. 
(arcsec) 

PA 

(°) 

2011 Mar. 26 

HiCIAO 

H 

12 

ADI 

52.7 

2.479 ± 0.016 

327.94 ± 0.39 

2011 May 22 

HiCIAO 

H 

21 

ADI 

56.9 

2.483 ± 0.008 

327.45 ± 0.19 

2011 Aug. 12 

IRCS 

V 

7 

DI 

— 

2.481 ± 0.033 

326.84 ± 0.94 

2011 Aug. 15 

IRCS 

V 

19 

DI 

— 

2.448 ± 0.024 

325.82 ± 0.66 

2012 Feb. 28 

HiCIAO 

K s 

38 

ADI 

42.5 

2.483 ± 0.015 

326.46 ± 0.36 

2012 Apr. 12 

HiCIAO 

J 

40 

ADI 

66.1 

2.487 ± 0.008 

326.54 ± 0.18 

2012 May 25 

IRCS 

V 

68 

ADI 

129.0 

2.499 ± 0.026 

326.14 ± 0.61 


Note. — The dates, used instruments, observation bands, total detector integration times 
(DIT), observation modes, and field rotation angles (FRA) are summarized for each observation 
run. The measured projected separations and position angles (PAs) are also described. 


H - band follow-up observations were performed on May 
22, 2011. The new U-band observations had both a 
longer total exposure time and better AO performance 
than the March observations, giving a higher signal-to- 
noise ratio in the final images. We obtained the K s - and 
J-band data in February and April of 2012. The seeing 
in the J-band observation was among the best ever mea- 
sured during SEEDS observations, and the AO worked 
very well, but those in iWband observations were worse 
than J- and i7-band observations. 

To correct the optical distortions of HiCIAO (and 
IRCS), we observed the globular cluster M5 using both 
instruments and compared the results to distortion- 
corrected images taken using the Advanced Camera 
for Surveys on the Hubble Space Telescope; the accu- 
racy and calibration for plate scale, orientation angle, 
and distortion of ACS/WFC camera, whose data were 
used for our calibration, have been previously reported 
(van der Marel et al. 2007). The distortion-corrected 
plate scale was 9.500 db 0.005 mas/pixel for HiCIAO, 
and its north orientation was 0.35 ± 0.02 degrees. Un- 
certainties in the distortion correction are less than 4 mas 
for HiCIAO within the central 10" x 10" of the field-of- 
view (FoV). We have included all of these uncertainties 
in our positional measurements for the planet detected in 
both HiCIAO and IRCS images, though they have little 
effect on our results. Photometric calibrations were per- 
formed by observing the central star itself with neutral 
density (ND) filters with transmissions of 0.590, 0.0628, 
and 1.14% for J-, i7-, and i^ s -bands, respectively. 

3.2. IRCS Observations 

To obtain data at a longer wavelength, we observed 
GJ 504 in the L'-band using IRCS on August 12 and 
15, 2011. We did not use either ADI or an occulting 
mask during these long- wavelength observations. In ad- 
dition, we obtained each frame by dithering the target 
within the field-of-view, which allowed us to better sub- 
tract the high thermal background (see below). In order 
to confirm whether we can obtain advantages for L'-band 
observational results by employing ADI, we performed 
additional L'-band observations in May 25, 2012. The 
longer total exposure time, compared to the 2011 ob- 
servations, as well as the usage of ADI, led to a higher 
signal-to-noise ratio. 

As in the case of HiCIAO, we observed M5 to obtain 
the calibration data for correcting the distortion of IRCS. 


After the distortion-correction, IRCS has the plate scale 
of 20.54 ± 0.03 mas/pixel, and its north orientation was 
tilted by 0.28 =b 0.09 degrees. Uncertainties in the dis- 
tortion correction are less than 8 mas for IRCS within 
the central 10" x 10" of the FoV. Because the L'-band 
magnitude of GJ 504 has not been measured, we cal- 
ibrated it using a standard star observed right before 
GJ 504 observations. As with HiCIAO observations, we 
then measured the flux of GJ 504 with an ND filter with 
a transmission of 0.67%. 

3.3. High- Contrast Data Reduction 
3.3.1. HiCIAO Data Reduction 

For the HiCIAO data reduction, we first removed stripe 
patterns (see Suzuki et al. 2010) emerged on each data 
frame, and subsequently performed the flat-fielding. Af- 
ter correcting the image distortion, we registered each 
dataset to the centroid of GJ 504 A. For the HiCIAO ob- 
servations, which used an opaque O'/ 4 diameter occulting 
mask, we estimated the relative frame centers by cal- 
culating cross-correlations of the frames as a function of 
positional offset. The A0188 and ADC kept the PSF ex- 
tremely stable and well- centered. For confirmation, we fit 
linear and quadratic functions to the frame-to-frame po- 
sitional offsets and estimated the positional drifts that re- 
mained after AO/ ADC correction. The systematic drifts 
were less than 10 mas for all observations except K s - 
band (17 mas), with residual rms scatters of ~2-10 mas. 
These random positional jitters can be attributed to our 
registration inaccuracies. 

We calibrated the absolute centers using a sequence of 
unsaturated, unmasked exposures taken with a ND fil- 
ter. For the i7-band data, the positional reference was 
acquired right after the end of the masked sequence. In 
the cases of the J- and iCg-band observations, we ob- 
tained the positional calibration data in a similar way 
as for the i7-band observations, but with the calibration 
observations interspersed through the science sequence. 
The positional errors for the unsaturated data ranged 
from ~ 2 to ^5 mas, where the lower quality PSFs for 
the i7-band data obtained in 2011 March observations 
or the Ag-band data in 2012 February caused the largest 
errors. The unsaturated data used a combination of ND 
filters and field lens that were different from the main 
science data, possibly introducing a small systematic off- 
set in position. To measure offsets introduced by the 
insertion of ND filters, we obtained images of a pinhole 
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mask at the focal plane of HiCIAO, with and without 
ND filters, and measured the offsets. To determine off- 
sets introduced by the use of a different field lens, we 
examined stars in images of the M5 globular cluster, ob- 
served both with and without the mask field lens. We 
measured these offsets to be small, < ^3 mas. The val- 
ues of the offsets introduced by the mask field lens were 
confirmed with three binaries in the SEEDS samples. We 
did not attempt to remove the offsets, but we have in- 
stead considered both of the offsets in the measurements 
of positions by adding systematic errors of 3 mas to the 
image registration uncertainties. In total, for iL-band ob- 
servations, image registration conservatively contributed 
~12 mas to our reported positional uncertainties for GJ 
504 b in March and mas in May. Those uncertainties 
are ~10 and ^6 mas in February (K s ) and April (J) of 
2012, respectively. 

We used the Locally Optimized Combination of Images 
(LOCI) algorithm (Lafreniere et al. 2007) to reduce the 
HiCIAO data. LOCI constructs a model PSF for each 
frame using the other frames in that dataset as refer- 
ences. The resulting PSF is locally a linear combination 
of the reference images, with the coefficients determined 
using a least-squares minimization. To avoid subtracting 
astrophysical sources, LOCI only uses reference frames 
with a minimum level of field rotation. For our reduc- 
tions of GJ 504, we required at least 1.2 ( =N a ) x FWHM 
/ R (x 180 / 7r deg) of field rotation, where the FWHM of 
PSFs on HiCIAO images are ^50-60 mas and R is the 
radial distance from GJ 504. As a result, any physical 
source was displaced by at least 1.2 FWHM between the 
working frame and the reference frames. We investigated 
if the signal-to-noise ratio was improved by applying N a 
= 1.5 to the reductions but the signal-to-noise ratios re- 
mained almost the same or were slightly reduced com- 
pared with the cases for N a = 1.2, and we thus carry out 
all analyses based on the images reduced with N a = 1.2. 

For the iL-band observations, our analysis of the data 
from March 26 detected a faint companion at a signal- 
to-noise ratio (SNR) of 5, where we calculated the noise 
as the standard deviation of the intensity in concentric 
annuli about GJ 504. Because of the longer integration 
time and better observing conditions for the observations 
in May 2011, we detected GJ 504 b with SNR = 9.0. We 
detected no other companion candidates (SNR < 5) in 
either observation. GJ 504 b is clearly confirmed with 
SNR of 18.4 at J-band, while the SNR of the iG-band 
observations is 4.6. 

The J, iL, and iG-band images for GJ 504 b are shown 
in Fig. 4. Also, the high-contrast two-color composite 
images is shown in Fig. 5, in which GJ 504 b is clearly 
visible 2'/48 northwest of GJ 504. 

3.3.2. IRCS Data Reduction 

For the IRCS data, we first corrected flat-fielding. Sub- 
sequently, we subtracted the thermal background and 
then used an algorithm (Galicher et al. 2011) similar to 
LOCI in order to further reduce the background. We es- 
timated the background value at each pixel and in each 
frame by dithering the images and using a linear com- 
bination of the frames in which the central source, GJ 
504, is sufficiently displaced. As in LOCI, we computed 
the coefficients for this linear combination using a least- 
squares fit. Eventually, we subtracted the resulting esti- 


mate of the local background whose variation is a func- 
tion of time. As in HiCIAO reductions, we subsequently 
corrected the distortion of IRCS images. 

The IRCS L'-band data were taken without a mask, 
and were saturated out to a radius of about 3-4 pixels 
(^60-80 mas, 0.6-0. 8 FWHM, where FWHM is the full 
width at half maximum of the PSF). We first performed 
relative registration by cross-correlating these saturated 
images with one another. We then centroided a set of 
unsaturated reference images, and cross-correlated the 
mean of these unsaturated PSFs with each saturated 
frame to measure the positional offset. Finally, we ap- 
plied this average offset to each of the saturated frames. 
The rms scatter of the offsets was ^6 mas for both obser- 
vations in 2011 August, and ^4 mas for 2012 May. We 
adopt these scatters as the registration uncertainties in 
L'-band data. 

We then combined these registered frames. We de- 
tected GJ 504 b with a signal-to-noise ratio of 3 and 4.3 
in the August 12 and August 15 data, respectively. Be- 
cause the data for 2012 May 25 were obtained with ADI, 
we could use LOCI for the PSF subtraction before stack- 
ing the registered data: the LOCI reduction parameters 
were set to be the same as the case of HiCIAO. In the re- 
sulting image, the companion was clearly detected with 
a signal-to-noise ratio of 8.0 (see Fig. 4). 

3.4. Positional Measurements of GJ 504 b 

We measured the position of GJ 504 b by fitting 2D el- 
liptical Gaussian functions to the final reduced data. Un- 
certainties in these Gaussian centroids range from 2 to 44 
mas, and dominate the total uncertainties in the position 
of GJ 504 b. Sometimes, LOCI introduces artificial shifts 
of the companion positions (Soummer et al. 2011). Mea- 
suring positions of artificial companions, which have the 
similar brightness and locations with the GJ 504 b, we 
have investigated how large offsets are added to the posi- 
tions of GJ 504 b. Because our adaptation for LOCI pa- 
rameters is not aggressive and the planet’s separation is 
sufficiently large from the star (r = 2 / /48), the astromet- 
ric biases are not large 6 . The positional measurements 
are converted to the projected separations and position 
angles (PAs) relative to the central star, and all of those 
are summarized in Table 2. Then, the uncertainties of 
measured separations and position angles were calculated 
by including the possible individual error sources of the 
image registration, distortion correction, plate-scale and 
orientation angle of each camera, and PSF-fit for GJ 504 
b. 

3.5. Photometry 

We performed aperture photometry on all of our ob- 
servations of GJ 504 b. For HiCIAO observations 
of the star GJ 504, we adopted magnitudes of 4.11, 
3.86, and 3.81 for its J-, iL-, and iG-band brightness 
(Kidger & Martfn-Luis 2003). The L'-band brightness 

6 Typically, averages of amplitudes for offsets are smaller than 3 
mas for all epochs of data, and the rms scatters are also sufficiently 
small compared to our assigned total positional errors. We do not 
add the possible systematic errors due to LOCI processing to our 
astrometric total errors since our measurements for the systematic 
offsets are degenerate to other error sources such as speckles or 
photon noise. However, if we include those in final astrometric 
errors, they have no impact on our conclusions. 
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Fig. 4. — J (~1.2 n m), H (~1.6 fi m), K s 2.2 /im), and L'-band 3.8 fi m) images of the newly discovered exoplanet, GJ 504 b. 
The four top panels show images reduced with the LOCI pipeline (a: J, b: H , c: K s , d: U). The corresponding signal-to- noise maps are 
shown on the four bottom panels (e: J, f: H , g: K s , h: Z/), in which the planet is detected with signal-to-noise ratios of 18.4, 9.0, 4.6, 
and 8.0, respectively. All signal-to-noise maps are shown at a stretch of [—5, 5]. In all panels, the star is located approximately in the lower 
left corner. North is up and East is left. 
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Fig. 5. — Discovery images of the exoplanet GJ 504 b. The central 6" X 6 ,x of two different high-contrast images were overlaid for this 
false color-composite image: orange represents iZ-band (~1.6 /i,m; Subaru/HiCIAO; May 22, 2011), and blue represents J-band (~1.2 /i,m; 
Subaru/HiCIAO; April 12, 2012). The J-band image was rotated by 0.9 degrees to compensate for the planet’s observed orbital motion 
(cf. Fig. 7). Panel a shows the intensity after suppressing flux from the central star. Panel b shows the associated signal-to-noise ratio 
(~9 in iJ; ~18 in J). The orbital radius of Neptune (~30 AU) is shown for comparison with our own solar system. The planet GJ 504 b 
is clearly visible as a white spot at a projected distance of 43.5 AU from the star GJ 504. The white color implies that the planet signal is 
persistent in both observations, setting it apart from the uncorrelated residual noise in each of the constituent images. 






Direct Imaging of GJ 504 b 


11 


Projected Separation (AU) 


o 

LO 


in 

S 

c 

o 

u 



Projected Separation (arcsec) 

Fig. 6. — Detection limits for J, H , K s , and L'-bands. Each panel shows 5cr detection limits (blue curves) estimated with reduced J, H , 
K s , and Z/-band data. The left vertical axes indicate measured contrasts relative to the primary at each separation. The corresponding 
detectable mass limits, at the system age of 300 Myr using COND models (Baraffe et al. 2003), are shown along the right axes. The bottom 
horizontal axes show angular separations projected from the primary in arcseconds. The top axis gives separations in AU. The contrasts 
for GJ 504 b are indicated by black squares in each panel, a, J-band detection under exceptionally good weather and seeing conditions; b, 
iJ-band; c, AT s -band; d, L'-band. 
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of GJ 504 was calibrated with a standard star, so that 
we obtained U = 3.94 db 0.09 mag for GJ 504. For all 
photometries of both the primary and the planet, we set 
the aperture radius to be 1.0 PSF FWHM, which ranged 
from 47 mas in J to 100 mas in L ' . 

By subtracting a (local) linear combination of frames 
from each image, LOCI can artificially reduce the flux 
of a point source (Lafreniere et al. 2007). We limit this 
effect by choosing a relatively large N a and estimate its 
magnitude using artificial point sources. We inject point 
sources in the HiCIAO data with projected separations 
and fluxes similar to those of GJ 504 b, measuring their 
magnitudes before and after applying LOCI. We mea- 
sured flux losses ranging from 13% to 23% for all LOCI 
reductions, where the typical scatter for the measure- 
ments of flux losses was 5%. We added this value to 
the final photometric errors. To verify this result, we re- 
reduced the HiCIAO data by subtracting a median PSF 
(Marois et al. 2006), but not applying LOCI. The result- 
ing fluxes for GJ 504 b had larger errors, but were consis- 
tent with their values in the final, LOCI-processed data. 
We estimated all of our photometric errors for the com- 
panion by calculating the rms scatter of aperture fluxes in 
blank-sky regions at the separation of GJ 504 b. For the 
iL-band data, these photometric errors were larger than 
the flux difference between the two sets of observations. 
For iL-band observations, the unsaturated frames used to 
calibrate the companion fluxes were obtained before and 
after the main science sequences with the mask, while the 
unsaturated integrations for J-, iL s -, and L'-band were 
conducted in intervals among the main sequences as well. 
The total number of unsaturated frames that were used 
for photometric calibration was more than 31 for each 
epoch observation, except in March 2011. We measured 
standard deviations of photometry of the primary using 
the unsaturated data; the worst value for the standard 
deviation was 22% in March 2011, while these values 
were less than 12% in all other observations. We include 
these variations in the final photometric uncertainties, 
in addition to the errors of photometry for companion 
and self-subtraction estimates. For L'-photometry, we 
had three sets of images, but did not combine the re- 
sult of 2012 with those derived from the data obtained 
in 2011, because the signal-to- noise ratios of the L' im- 
ages in 2011 were significantly worse than for the 2012 
images; in the former case, no ADI was used and the 
observing conditions were very poor, with airmasses as 
high as 3.5 at maximum. Therefore, our adopted value 
for the L' photometry is based exclusively on the data 
obtained in 2012. The photometric measurements of GJ 
504 b at each band are summarized in Table 1, with the 
iL-band flux computed using the weighted mean of the 
March and May values. Furthermore, the detection lim- 
its estimated from the reduced data are shown in Fig. 
6 . 

4. RESULTS 

4.1. Astrometry 

To confirm that GJ 504 b is gravitationally bound to 
the host star, we tracked its proper motion over the 

course of more than a year, from March 2011 through 
May 2012 (Fig. 7). Fig. 7 allows to clearly visually 
distinguish the measured positions from the track ex- 


a Projected E/W separation (") 

0.8 0.9 1.0 1.1 1.2 1.3 1.4 



Fig. 7. — Astrometric analysis of GJ 504 b. The measured po- 
sitions of the planet relative to its parent star at the seven epochs 
of observation are shown as plus signs (size corresponds to error 
bars). A background star would follow the dotted trajectory, with 
star symbols marking the position at the epochs of observation. 
More information on the orbital simulation can be found in the 
Appendix A. Diamond symbols mark the position of the planet 
along the most likely orbit at the epochs of observation. Plus, star, 
and diamond symbols are color-coded for epoch, a, View of the 
planet’s motion compared with the track expected if GJ 504 b were 
a background star; b, zoomed-in view of the observed positions of 
GJ 504 b. 

pected for a background object. However, we also car- 
ried out a x 2 test to derive a statistical confidence that 
the source is a background star. This yielded a value 
of x 2 = 1,340 (12 degrees of freedom). Therefore, the 
X 2 * * test on the astrometry decisively excludes GJ 504 b 
as an unrelated background source. On the other hand, 
if we assume GJ 504 b to be an interloping foreground 
object, its color ( H — L' = 3.3 mag) would rule out a 
star anyway, leaving only the possibility of a cool brown 
dwarf (T < 600 K) within a distance of ^30 pc. Even un- 
der optimistic assumptions about the brown dwarf mass 
function (Kirkpatrick et al. 2011), less than one such ob- 
ject exists per 10 square degrees, and the probability of 
a chance alignment within 3" of GJ 504 is no more than 

io- 7 . 

Furthermore, in order to demonstrate that the mea- 
sured motion of GJ 504 b relative to its parent star GJ 
504 is indeed consistent with the expected orbital mo- 
tion of a gravitationally bound planet, we have fitted a 
Keplerian curve to its measured positions relative to the 
star, and constrained its orbit (see Appendix A). Here 
we assumed that the star has 1.22 M 0 . We show the 
best-fit curve in Fig. 7. The best-fit y 2 was derived to 
be 11.7 (8 degrees of freedom), which confirms that the 
measured positions relative to the central star are con- 
sistent with dynamical orbital motion around the star 7 . 
These analyses robustly prove that GJ 504 and GJ 504 
b comprise a physically associated system. 

4.2. Mass Estimate 

Using the luminosity and age of GJ 504 b, we thus 
proceed to estimate its mass based on the hot-start model 
(Baraffe et al. 2003) (see Fig. 8). We then calculate 

7 The position measured in Aug. 15 of 2011 deviates from the 
best-fit orbit, which increases the % 2 - That measurement was ob- 
tained without ADI at high airmass. 
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Fig. 8. — The mass of GJ 504 b, as derived from the theoretical 
relation between the planetary age and luminosities for various 
wavelengths. Estimated age (Myr) and near-infrared luminosities 
(absolute magnitudes) for GJ 504 b are plotted over cooling curves 
for planets of various masses calculated from the hot-start model 
(Baraffe et al. 2003); a, L'-band (~3.8 pm) luminosity; b, K s - 
band (~2.1 pm); c, if-band (~1.6 pm); d, J-band (~1.2 pm). 
The vertical size of each box shows the estimated luminosity range 
(~ lcr) in each band. The gyrochronological age (indicated by the 
gray box) is the most direct age estimate among all age indicators 
available for GJ 504; it implies a planet mass of 4 ^'q Mj up . The 
open box with the solid black line boundary indicates the wider age 
range encompassing all of the adopted age estimates, and implies 
a mass of 4 +-['q Mj up . 


a mass by comparing luminosity of GJ 504 b at each 
band to the models, and eventually average the mass 
estimates at J, iL, K s , and L'-bands, which provides a 
mass of the planet, and its error range, which includes 
the errors from age and photometry. Any bias in the 
mass estimate caused by focusing on the photometry in 
a particular band should be mitigated when taking the 
average of the photometry of all bands. For the upper 
bound, we estimate the mass of GJ 504 b from the la 
bright end of the photometry at each band, using hot- 
start models for an age of 510 Myr. The collective upper 
limit is then acquired by averaging the estimates at each 
band. The lower limit is calculated in an equivalent way, 
using the lcr faint end photometry and an age of 100 Myr. 
Likewise, the best estimate uses the mean photometry 
in each band and an age of 160 Myr. As a result of 
this procedure, a mass estimate of 4 ^'q Mj up is derived. 
The age uncertainty dominates the errors. If we adopt 
the gyrochronological age of 160^gQ Myr, the most direct 
of our age indicators, the mass of GJ 504 b is constrained 
to lie between 3 and 5.5 Mj up . This is well below 13.6 
Mj up , the deuterium-burning threshold commonly used 
to divide planets and brown dwarfs (Marois et al. 2008; 
Burrows et al. 1997) and among the lowest masses for 
exoplanets discovered by direct imaging. 

Adopting the iG-band luminosity may result in a bias 
in the mass estimation, since this wavelength range is 
particularly sensitive to atmospheric parameters such 
as surface gravity and metallicity (Fortney et al. 2008; 
Spiegel & Burrows 2012; Janson et al. 2011). Indeed, 
Fortney et al. (2008) predict that raising the metallicity 
to 5 times Solar metallicity results in a iG-band bright- 
ening of a factor 2-6 for a planet with temperature and 
surface gravity values in the range of what is expected 


for GJ 504 b. Since the star GJ 504 has a metallicity 
above the Solar value, it is very likely that the planet’s 
metallicity is also enhanced. For reference, if we neglect 
the iG-band photometry, the best-fit mass of GJ 504 b 
becomes 3.5 Mj up , and is confined to be less than 8 Mj up 
based on the adopted age-range. 

The robustness of mass estimates for directly imaged 
exoplanets is usually limited by theoretical models used 
to derive age- luminosity relations. One uncertainty in 
these relations is the choice of initial conditions, which 
influences the estimation of a planetary mass. Since 
GJ 504 is older than 100 Myr, its inferred mass is less 
sensitive to the uncertainty in initial conditions for the 
employed model (Mar ley et al. 2007; Spiegel & Burrows 

2012) than the previously imaged planets which have 
younger age than GJ 504 b. To test the sensitivity of the 
G J 504 b mass to uncertainties in the hot-start model, we 
compare the near-infrared luminosities of GJ 504 b with 
a suite of models with ages between 100 and 510 Myr, 
initial thermal states ranging from hot to extremely cold 
(initial specific entropies from 8.0 to 12-13 &b baryon -1 ), 
masses between 1 and 20 Mj up , varying cloud properties, 
and metallicities between one and three times the Solar 
value (Spiegel & Burrows 2012). Even when adopting 
an age of 100 Myr, for which the largest differences arise 
between the hottest- and coldest-start models, the dif- 
ference in the mass estimates between the two cases is 
typically ^4 Mj up , which is calculated as the median 
of all estimates for each different model atmosphere and 
each band (see Fig. 9). At older ages, the difference de- 
creases even further. By contrast, if we use HR 8799 b 
(20-50 Myr; Marois et al. 2008, 2010) as an example of 
planets younger than 50 Myr, the corresponding differ- 
ence is higher than 10 Mj up , which is significantly larger 
than the case of GJ 504 b. Thus, in quantitative terms, 
the dependency on initial conditions for GJ 504 b is in- 
deed rather week. This also allows to confine the mass of 
GJ 504 b to within the planetary mass regime, essentially 
regardless of initial conditions (See Fig. 10). The mass 
inferred from each of the J-, i7-, and L'-bands individ- 
ually is practically always below the deuterium burning 
limit of ^14 Mj up , apart from the most extreme cases in 
which both the coldest initial conditions and ages near 
the oldest end of the age range are adopted simultane- 
ously. Note that these extreme cold- start models may be 
unrealistic for giant planets (Spiegel & Burrows 2012). 

5. DISCUSSION 

Here, we discuss the unique properties of GJ 504 b 
through comparison with other imaged planets. Further- 
more, the possible origin of GJ 504 is discussed. 

5.1. Previously Imaged Planet Candidates 

Some of the previously imaged exoplanets were dis- 
covered around young (10-50 Myr), relatively massive 
stars like HR 8799 (1.5 Ms un ; 20-50 Myr; Marois et al. 
2008, 2010) and /? Pictoris (1.75 Ms un ; 10-20 Myr; 
Lagrange et al. 2010), with masses estimated at 5-9 
Mj up , assuming hot-start conditions. Recently, SEEDS 
detected a ^13 Mj up companion around the B9-type 
star Kappa And (2.5 Ms un ; 20-50 Myr; Carson et al. 

2013) . An exoplanet candidate around another mas- 
sive star, Fomalhaut (2.0 Ms un ; ^200 Myr; Kalas et al. 
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Mass (M Jup ) 

Fig. 9. — Differences between mass estimates from cold-start models and hot-start models for 100 Myr old planets. The mass-luminosity 
relations modeled in Spiegel & Burrows (2012) are illustrated with dashed and solid lines for hot-start (entropy = 13 &b baryon -1 ) and 
cold-start (8 &b baryon -1 ) initial conditions, respectively. We show the relations for four atmosphere types on each panel (a: J, b: if, c: 
K s , d: L'): “hybrid” means patchy cloud model and 1 or 3 S means one or three times solar-abundance model. The gray regions correspond 
to the absolute magnitudes of GJ 504 b in each band and their error ranges (Tier). The lower age limit (100 Myr) of GJ 504 is chosen in 
this illustration since that is the stage at which the largest differences occur between hot-start and cold-start conditions. As an example, 
the arrow in panel b indicates the difference in mass arising from different initial conditions when the IS cloud- free model atmosphere is 
adopted. The difference is ~4 Mj up , which is also the median result for the full set of atmospheres and photometric bands. 


2008), has been detected at optical wavelengths. The 
images were reanalyzed by Galicher et al. (2013) and 
Currie et al. (2012), who recovered the point source at 
an additional wavelength. In addition, new optical im- 
ages were recently presented with a report of an eccen- 
tric orbit (Kalas et al. 2013). Very recently, another 
imaged planet candidate has been reported around the 
young star HD 95086 (10-17 Myr; ^1.6 Ms un ), but a 
robust second epoch confirmation still needs to be ac- 
quired (Rameau et al. 2013). Other giant planet candi- 
dates have also been imaged at wide separations around 
very young (<10 Myr) stars (See e.g., Kuzuhara et al. 
2011). Among these, 1RXS J1609 b (Lafreniere et al. 
2008) orbits a pre-main sequence solar-mass star at a 
wide separation of 330 AU. 

For all the aforementioned companions, their masses 
have been inferred from comparisons of their infrared lu- 
minosity with hot-start models. However, due to the 
youth of those systems, they are highly affected by hot- 
versus cold- start conditions. The cold- start models imply 
that even HR 8799 b, which may have the lowest masses 
among them, would be more massive than ^14 Mj up 
(See Fig. 10): its mass is estimated to be ^20 Mj up or 
higher (except at J-band). Therefore, under cold-start 
assumptions, all the exoplanets around these young stars 
have much higher masses than currently inferred, and are 


pushed into the brown dwarf mass regime (>~14 Mj up ). 
One special case is the candidate around Fomalhaut, for 
which a mass below 3 Mj up has been reported, but deep 
infrared non-detections indicate that the detected light 
from the candidate cannot arise from a planetary sur- 
face (Janson et al. 2012); hence, its physical nature re- 
mains unclear. Regardless of the caution for the hot-start 
models, a few independent constraints, such as dynami- 
cal stability analysis in the case of the HR 8799 system 
(e.g., Fabrycky & Murray-Clay 2010; Marois et al. 2010; 
Sudol & Haghighipour 2012), suggest that the hot-start 
models may be realistic. 

As discussed above, the mass estimates for directly- 
imaged planets suffer from uncertainties of cooling mod- 
els for giant planets until future calibrations may be ob- 
tained by comparing model-dependent masses with im- 
proved dynamical mass estimates. However, due to its 
old age, GJ 504 b is less dependent on the models’ uncer- 
tainty related to the hot versus cold conditions, compared 
with previously imaged planets. 

5.2. Color of GJ 504 & 

We can compare the magnitudes and colors of GJ 
504 b to similar objects in the literature in order 
to get a sense of the physical properties of this 
planet. For this purpose, it has been placed on color- 
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Fig. 10. — Comparison of multiple wavelength magnitudes of GJ 504 b with cold-start model (Spiegel & Burrows 2012) predictions at the 
planet/brown dwarf boundary. In each panel, the cooling tracks predicted by the coldest-start models for a 14 Mj up planet are plotted in 
the age range of 20-500 Myr. Different panels correspond to the different bands (a: J, b: i7, c: K s , d: L'). In addition, we show the cooling 
curves for four atmosphere types, where “hybrid” means patchy cloud model and 1 or 3 S means one or three times solar-abundance model 
(Spiegel & Burrows 2012). The age range of GJ 504 b is indicated by solid and dotted lines. The solid lines and dotted lines correspond 
to the age range inferred with gyrochronology (160^gQ Myr) or chromospheric activity (330 d= 180 Myr), respectively. For comparison, we 
plot the brightness and age of HR 8799 b (30^q Myr) (Marois et al. 2008; Metchev et al. 2009) in each panel. 
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Fig. 11. — Color-magnitude diagrams for GJ 504 b. Shown for comparison are LT-type field brown dwarfs as black dots (Leggett et al. 
2010), and a sample of planet and brown dwarf companions from the literature. Panel a shows J-band luminosity (J-band absolute 
magnitude) versus J — H and panel b shows if-band luminosity versus H — K s ( H — K c for GJ 758 B). Directly imaged planets are each 
labeled. 
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magnitude diagrams (Fig. 11) along with LT-type field 
brown dwarfs (Leggett et al. 2010) and HR 8799 bcde 
(Marois et al. 2010; Metchev et al. 2009; Skemer et al. 
2012), Kappa And b (Carson et al. 2013), 1RXS J1609 b 
(Lafreniere et al. 2008), 2M 1207 b (Chauvin et al. 2004; 
Mohanty et al. 2007; Skemer et al. 2011), and GJ 758 B 
(Thalmann et al. 2009; Janson et al. 2011). 

Given the location of GJ 504 b in the color-magnitude 
diagrams, it is clear that it is much colder than previous 
planet detections. Indeed, an application of JHK S L' lu- 
minosities of GJ 504 b to those predicted by the cooling 
models of Baraffe et al. (2003) implies an effective tem- 
perature of 510^20 K for the planet, where we estimated 
the best estimate and lower/upper limits in the same way 
as for the mass estimate of GJ 504 b (see section 4.2). 
The inferred effective temperature is significantly colder 
than the values of 800-1800 K inferred for previously 
imaged planets (Marois et al. 2008; Chauvin et al. 2004; 
Bonnefoy et al. 2011). The extrasolar giant planets other 
than GJ 504 are primarily located around the so-called 
LT transition, which corresponds to the transition from 
cloudy (L-type) to clear (T-type) atmospheres. GJ 504 
b is consistent with having a largely clear atmosphere 
analogous to the bottom of the T-type sequence, or the 
top of the yet colder Y-type sequence. It deviates from 
this sequence in H — K s , due to over- luminosity in K s - 
band, which is probably due to a low surface gravity or 
enrichment in heavy elements, both of which can repro- 
duce such an effect (Burrows et al. 2006). These factors 
have been shown to possibly explain the same trend seen 
in GJ 758 B (Janson et al. 2011), which is an older and 
more massive companion to a Sun-like star. We note that 
GJ 504 b was observed with the iG-band filter and GJ 
758 B with the iG-band filter, which prevents a one-to- 
one comparison. However, models incorporating heavy 
element enrichment appear to be a promising option for 
reproducing the full set of magnitudes of GJ 504 b. 

It has been proposed that planetary candidates with 
hot atmospheres (T e ff ~ 1000 K) around HR 8799 
or 2M 1207 may have thick clouds (Barman et al. 
2011a; Madhusudhan et al. 2011; Skemer et al. 2011; 
Mar ley et al. 2012) although their cloud properties are 
under discussion (Barman et al. 201 la, b). By contrast, 
the relatively blue J — H (= -0.23) color of GJ 504 b is 
consistent with T-type brown dwarfs in the same temper- 
ature range (Leggett et al. 2010), which are representa- 
tive of less cloudy atmospheres that could occur naturally 
in this temperature range (Allard et al. 2001). In addi- 
tion, models examining cloud clearing as a function of 
temperature and surface gravity in Mar ley et al. (2012) 
do predict that a temperature of ^500 K and log# of 
~4, as expected for GJ 504 b (see Table 1), should lead 
to a clear atmosphere, while the properties of previously 
imaged planets place them in a cloudy regime. The cold 
and perhaps less cloudy atmosphere of GJ 504 b thus 
places it in a physically distinct state from the hotter 
and cloudier atmospheres of previously imaged planets, 
and should be highly interesting for further atmospheric 
studies in the future. These properties imply that GJ 
504 b will become a benchmark object for the study of 
exoplanet atmospheres. 

5.3. Implications for the Origin of GJ 504 b 


GJ 504 b is a giant planet that lies at a projected 
separation of 43.5 AU. Two major scenarios could ex- 
plain its formation: the core-accretion (CA) model 

(Pollack et al. 1996) and the gravitational instability 
(GI) model (Durisen et al. 2007). Here, we discuss 
whether the CA or GI can account for the origin of GJ 
504 b. 

5.3.1. Core Accretion 

In the CA model, the core of a giant planet accretes 
planetesimals and grows to a critical mass (~10 Earth 
masses), at which point the core begins to rapidly ac- 
crete gas directly from the protoplanetary disk. It is 
difficult for the CA model to explain the formation of 
giant planets in situ beyond 30 AU. A recent study 
(Dodson-Robinson et al. 2009) suggested that even a lo- 
cally stable (but globally unstable) disk with an order 
of magnitude higher mass than the minimum mass so- 
lar nebula (Hayashi 1981) around a G-type star like GJ 
504 was unable, in simulations, to form giant planets 
beyond 30 AU within the lifetime of the protoplanetary 
disk (Dodson-Robinson et al. 2009). Therefore, the pres- 
ence of GJ 504 b, as well as other directly-imaged planets 
such as HR 8799 bcde, remains a particular challenge for 
the conventional core accretion theory. In this implica- 
tion, we should note such numerical investigations as con- 
duced by Dodson-Robinson et al. (2009) have typically 
neglected some theoretical processes that may enable the 
cores of giant planets to grow in the shorter timescale: 
effective damping of small planetesimal’s eccentricity due 
to gas drag (Rafikov 2004) or efficient capture of plan- 
etesimals due to atmospheric drag (Inaba & Ikoma 2003) 
may accelerate the accumulations of planetesimals in 
such the outer region, so that the in situ formation of 
GJ 504 b may become possible (Rafikov 2011). How- 
ever, even the models incorporating these processes may 
not necessarily be promising for producing a planet more 
massive than the critical mass beyond 30 AU at least for 
disks that are not particularly massive, as demonstrated 
by the simulations of Kobayashi et al. (2011). 

Even if the in situ formation of a giant planet is 
difficult in such outer regions, models incorporating 
planet scattering can account for the origin of GJ 504 
b. If multiple planets are formed by core accretion in 
the inner regions, one or more of them could be scat- 
tered outward due to dynamical interactions among the 
planets (Nagasawa et al. 2008) or between the planets 
and the protoplanetary disk (Paardekooper et al. 2010; 
Crida et al. 2009). GJ 504 b may have been pushed out- 
wards by as-yet undiscovered companions. Multiple gi- 
ant planets may be common in metal-rich protoplanetary 
disks (Fischer & Valenti 2005), because the high metal- 
licity may enhance solid materials and planetesimals in 
their protoplanet ary disks, resulting in rapid growth of 
planetary cores. Eventually, outward migration may fre- 
quently occur in the metal-rich systems. Hence, the rel- 
atively high metallicity of GJ 504 system (see Table. 1) 
may be a preferential property to explain the origin of 
GJ 504 b under the outward migration hypothesis. 

5.3.2. Gravitational Instability 

In contrast to CA, the GI model can form massive 
planets in situ at large radii. In this model, a massive 
protoplanetary disk becomes gravitationally unstable in 
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its outer regions (> ^50 AU) as the cooling time be- 
comes short relative to the local dynamical timescale 
(Durisen et al. 2007). The outer disk then fragments, 
collapsing directly into one or more giant planets. The 
instability can also occur between 20 and 50 AU, with 
continuous inflows of gas into the magnetically inactive 
region inducing fragmentation on a timescale of 10 3 yr 
(Machida et al. 2011). However, the planet would have 
to escape the very rapid inward migration expected in 
a turbulent disk (Baruteau et al. 2011). In addition, if 
multiple massive planets are formed in the GI disk, they 
would still be scattered and eventually ejected from the 
systems, as in the case of CA. Therefore, GJ 504 b would 
be a planet surviving against falling onto the central star 
or ejections form the system, if it was formed through 
the GI process. We should note that fragmentation could 
also be suppressed in a metal-rich disk like that expected 
for GJ 504; such a disk would cool less efficiently due to 
the higher opacities (Cai et al. 2006). Also, the signifi- 
cant gas accretion onto a clump formed according to the 
GI process may prevent the planet from remaining less 
massive (Zhu et al. 2012; Boley et al. 2010), disallowing 
the occurrence of a low-mass planet like GJ 504 b. 

5.3.3. Summary for Discussions of the Origin of GJ 50 4 b 
and Future Prospect 

Based on the available data, we cannot conclusively 
determine which process led to the formation of GJ 504 b. 
However, it is interesting to note that GJ 504 b is the first 
wide-orbit giant planet detected around a demonstrably 
metal-rich star. As discussed above, the high metallicity 
may enhance massive planetary cores, while preventing 
the effective cooling of a protoplanetary disk, and thus 
be more hospitable for CA formation than GI formation. 
Thus, given these measurable properties, planets such as 
GJ 504 b may provide important clues for understanding 
the formation and evolution of giant planets. 

For clarifying the origin of GJ 504 b, the further 
observations may be crucial. Intriguingly, the CA 
model suggests the presence of unseen massive planets 
at smaller angular separations (Nagasawa et al. 2008; 
Chatterjee et al. 2008), making GJ 504 an outstand- 
ing candidate for follow-up observations by current and 
future instruments. 7V-body simulations for the dy- 
namically unstable gas-free planetary system predict 
that such inner companions should orbit at a few AU 
from the host 8 and have masses equal or greater than 
the corresponding outer planets (Nagasawa et al. 2008; 
Chatterjee et al. 2008). Therefore, the inner counter- 
part for GJ 504 b would have to be more massive than 
^ 4 Mj up if it exists, although we note that radial ve- 
locity observations have found few planets in this mass 
range at semi- major axes of a few AU (Mayor et al. 2011) 
and the remaining gasses in the disk may affect the scat- 
tering of planets, changing the configurations of plan- 
ets predicted by the simulations for the gas-free cases 
(Moeckel & Armitage 2012). 

Also, further monitoring of the orbit may be impor- 
tant, since scattered planets preferentially have high ec- 
centricities (Chatterjee et al. 2008). Furthermore, mod- 
eling the atmosphere of GJ 504 b may provide another 

8 The current detection limit cannot explore giant planets in 
such the inner regions (see Fig. 6). 


piece of evidence to test whether it was produced by the 
CA process. The enhancement of heavy elements in the 
planetary atmosphere is a natural outcome of the CA 
process (Chabrier et al. 2007). It is possible to examine 
the abundances of such elements in the spectral energy 
distribution of GJ 504 b through the use of atmospheric 
models. For this purpose, multiple wavelength photo- 
metric observations as well as the spectroscopy can play 
an important role. 

Further observations of the GJ 504 system will con- 
strain the origin of GJ 504 b and may lead to a better 
understanding of the origin of planetary systems, includ- 
ing our Solar System. Future instruments with direct 
(e.g., GPI, SPHERE, SCExAO) or indirect (e.g., IRD; 
Tamura et al. 2012) techniques that can explore more 
inner region may be able to find a promising evidence to 
clarify the origin of GJ 504 b. 

6. CONCLUSION 

As part of the SEEDS direct imaging survey, we have 
detected a giant planet around the Sun-like star GJ 504. 
The star has a spectral type of GO and an approximate 
mass of 1.2 Ms un - Its age has been conservatively es- 
timated as 160lgo° Myr based on a combination of gy- 
rochronology and chromospheric activity. Gyrochronol- 
ogy alone, which is the most direct and likely the most 
reliable age estimator, places the age at the lower end of 
this range, at 160^0 Myr. Observations over one year 
baseline provided confirmation that the detected planet 
is orbiting the star. Compared with previously imaged 
exoplanets, GJ 504 b has a number of interesting fea- 
tures: 

1. The age of this planet is the oldest among all di- 
rectly imaged planets. 

2. Its mass is estimated to be Mj up . This is 

among the least massive of imaged exoplanets. Be- 
cause of its old age, the mass estimate is less depen- 
dent, compared to other imaged exoplanets, on the 
uncertainties of initial conditions in models used to 
convert its infrared luminosity into a mass. 

3. GJ 504 b is the first imaged giant planet on a wide 
orbit (> 5 AU) around a G-type main-sequence star 
with a mass near the solar mass. 

4. The planet’s projected separation from the pri- 
mary is about 43.5 AU, which is comparable to 
the widest-orbit planets around massive stars such 
as HR 8799. 

5. The effective temperature of the planet (510^0 K) 
is the coldest among all previously detected giant 
planets. 

6. The blue J — H color of -0.23 suggests that it has 
a less cloudy atmosphere than other imaged exo- 
planets. The cold temperature and color place the 
planet in a novel parameter space of atmospheres 
for exoplanets. 

7. GJ 504 b is the first giant planet on a wide orbit 
discovered around a demonstrably metal-rich star. 
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The performed observations alone cannot conclusively 
uncover the origin of GJ 504 b. Further observations 
of this system will enable more direct comparisons with 
our own Solar System, and will help unveil the formation 
history of giant planets in the outer disk. 
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APPENDIX 

MONTE CARLO ORBITAL ANALYSIS 

To demonstrate that the measured motion of GJ 504 b relative to its parent star GJ 504 is indeed consistent with the 
expected orbital motion of a gravitationally bound planet, we have run a Monte Carlo simulation following Janson et al. 
(2011). A large number of physically plausible orbits is randomly generated and compared to the astrometric data 
points. The maximum of the y 2 function is typically ill-defined for an astrometric dataset covering only a small fraction 
of an orbit: an entire family of orbits spanning a large range of semi-major axes match the data equally well. Thus, 
rather than ascribing an undue amount of significance to the simulated orbit with the smallest y 2 value, we consider 
the entire well- fitting family of orbits with y 2 < Xmin + 1. The plots in Fig. 12 summarize the properties of this 
well-fitting orbit family. 

To characterize this distribution numerically, we determine the weighted median and the weighted 68% range for 
each parameter of the orbital simulation (semi- major axis a, eccentricity e, inclination i, longitude of ascending node 
D, and argument of periastron cj). We calculate the statistical weight W of each orbit as the mean orbital velocity 
divided by the local orbital velocity at the epoch of observation. This represents the fact that an eccentric planet 
orbits slower, spends more time, and is therefore more likely to be observed at apastron than at periastron. 

However, some of the histograms are heavily skewed (e.g., log a) or even periodic (e.g., D); as a result, the median 
badly represents the true behavior of the entire well-fitting orbit family. For this reason, we further define the most 
likely orbit (MLO). 

For each simulated orbit n, we calculate as a measure of likelihood L n the product of the histogram values of the 
log a, e, i, D, uj histograms in the bin into which the orbit falls, as well as the statistical weight W: 


L n = histio ga (loga n ) • hist e (e n ) • histi (i n ) • hist o(D n ) • hist w (a; n ) • W. (Al) 

The MLO can then be defined as the orbit with the highest measure of likelihood, Lmlo = max n L n . 
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and the weighted 68% interval (dotted), whereas a blue dash-dotted line marks the most likely orbit. 
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